Strong-coupling expansions for the topologically inhomogeneous Bose-Hubbard model 



P. Buonsante/ V. Penna,^ and A. Vezzani^ 
^ Dipartimento di Fisica, Politecnico di Torino and I.N.F.M, 



o 
o 

(N 
> 

o 



o 



a 
o 

> 

o 
o 



I 

O 

o 



X 



Corso Duca degli Abruzzi 24 - 1-10129 Torino (ITALIA) 

^ Dipartimento di Fisica, Universitd degli Studi di Parma and I.N.F.M., 

Parco Area delle Scienze 7/a 1-43100 Parma (ITALIA) 

(Dated: February 2, 2008) 

We consider a Bose-Hubbard model with arbitrary hopping term and provide the boundary of 
the insulating phase thereof in terms of third-order strong coupling perturbative expansions for the 
ground state energy. In the general case two previously unreported terms occur, arising from tri- 
angular loops and hopping inhomogeneities, respectively. Quite interestingly the latter involves the 
entire spectrum of the hopping matrix rather than its maximal eigenpair, like the remaining per- 
turbative terms. We also show that hopping inhomogeneities produce a first order correction in the 
local density of bosons. Our results apply to ultracold bosons trapped in confining potentials with 
arbitrary topology, including the realistic case of optical superlattices with uneven hopping ampli- 
tudes. Significant examples are provided. Furthermore, our results can be extented to magnetically 
tuned transitions in Josephson junction arrays. 
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I. INTRODUCTION 

The Bose-Hubbard (BH) model describes a gas of in- 
teracting bosons hopping across the sites of a discrete 
structure. It was originally introduced^ to study the 
the superfluid-insulator quantum phase transition in liq- 
uid helium, and later argued to describe the physics of 
an array of coupled Bose-Einstein condensates (BECs)fSi 
Such proposal has been strikingly substantiated by a 
recent breakthrough experiment,^ where the superfluid- 
insulator transition was observed in a ^"^Rb BEC array. 

This result greatly reinvigorated the studies on 
the quantum phase transition characterizing the BH 
model4iS'SiS Indeed, the ongoing progress in conden- 
sate trapping techniques promises an unprecedented con- 
trol on the parameters of the experimental realizations of 
such model. A simple example in this respect is supplied 
by the quite direct relation between the intensity of the 
laser beam trapping the condensates at the sites of the 
so-called optical lattice and the tunneling interaction be- 
tween neighbouring condensates, accounted by the hop- 
ping amplitude in the relevant BH Hamiltonian>2i^ 

Recently, several proposals have been brought forward 
concerning the realization of BEC arrays characterized 
by non-trivial layouts. Suitable laser setups, possibly 
with different wavelengths, can be used to generate op- 
tical superlatticesAifliiiiiSiH Techniques borrowed from 
holographyi^ or from semiconductor technology""^^ have 
been proposed for the creation of optically or magneti- 
cally confined inhomogeneous BEC arrays. 

The BH model is also commonly used to describe 
Josephson junction arraysi^ (JJAs), which can be easily 
arranged according to no n- trivial layouts, including in- 
homogeneous onesii2ii£ In this case the phase transition 
can be observed measuring different samples or adjusting 
an external magnetic field applied to the same sample.i& 

In the light of this, the need arises for the study of a 



BH model characterized by a generic hopping term. A 
further point of interest in this respect is provided by 
the unexpected features induced by topological inhomo- 
geneities in the properties of discrete boson systemgiS 
even in the absence of boson interactioniS2i2i 

For the study of the zero-temperature phase dia- 
gram of the BH model several techniques have been 
adopted. Qualitatively correct results are provided 
by numericafeSSi and analyticalSiSiiS^ mean-field ap- 
proaches. As it is well known, even for relatively small 
systems the formidable size of the Hilbert space rules 
out the route of direct diagonalization in the study 
of the ground state properties of the BH Hamilto- 
nian, determining the features of the zero-temperature 
phase diagram. In the case of homogeneous systems, 
clever numerical methods have been devised that pro- 
vide quite satisfactory results. However, since they 
rely on features such as self-similarit>s2& or translational 
invariance^ their generalization to topologically inho- 
mogeneous structures is at least problematic. Of course, 
the superfluid-insulator transition on such structures 
can be studied exploiting the flexibility of Quantum 
Monte Carlo simulations, that provide the most gen- 
eral and reliable numerical approach. Indeed, this tech- 
nique was successfully applied to the case of regular 
lattice a^^i^^ , possibly in the presence of inhomogeneities 
arising from random22i2L2Si2^ or harmonic"''^ local po- 
tentials. Recently, some results addressing inhomoge- 
neous hopping have been reported»i^ As to analytical 
approaches, the strong coupling perturbative results re- 
ported in Ref. |3J allow a quite satisfactory description of 
the zero-temperature phase diagram of the BH model on 
regular bipartite lattices, possibly with disorder. 

Here , w e extend the third-order perturbative results 
of Ref. |3J to a generic structure, possibly non-bipartite 
and/or topologically inhomogeneous. We remark that 
commonly investigated regular lattices, such as the 



triangular^-'-- or the Kagome- ^ lattice, are non-bipartite. 
We show that in the general case two further significant 
third-order contributions appear in the perturbative ex- 
pansion of the ground state energy of the system, and of 
the boundaries of the Mott lobes thereof. Furthermore 
we show that terms arising from inhomogeneous hopping 
amplitudes appear in the expansion of a physically inter- 
esting quantity such as the local density of bosons already 
at the first perturbative order. We also supply some ex- 
amples where the new terms play a significant role. 



II. THE BOSE-HUBBARD MODEL ON A 
GENERIC STRUCTURE 

The Hamiltonian for the pure BH model on a generic 
structure comprising M sites has the form 
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where Qj (aj") is the boson annihilation (creation) oper- 
ator at site j, TLj = a^aj is the relevant number oper- 
ator, [/ > is the boson repulsion strength and t > 
is an overall scaling factor controlling the magnitude of 
the possibly inhomogeneous hopping matrix, t. Since the 
latter contains information about the connectivity of the 
underlying structure, it is ij^i = Q <^ Ajh = 0, where A is 
the so-called adjacency matrix, whose generic entry Ajh 
equals one only if sites j and h are adjacent, and zero 
otherwise. Hence topological inhomogeneity can come in 
the form of a non-trivial connectivity, i.e. a site depen- 
dent coordination number Zj — X^/i^j'n and/or in the 
form of uneven hopping amplitudes i^ . We remark that 
this applies also to periodic structures, such as the inho- 
mogeneous Kagome lattice in Fig.|5](C.) or the T^ lattice 
of Ref . llSl Note that Hamiltonian ^ commutes with the 
total number of bosons, N — X]i=i %• Hence the prop- 
erties of its ground state \^) — i.e. the state minimizing 
(^|7f — /iiVj?/;), where /i is the chemical potential control- 
ling the boson population — can be conveniently studied 
considering one number eigenspace at a time. 

The competition between the on-site repulsion and the 
kinetic energy — proportional to U and r, respectively — 
produces a well-known zero-temperature phase diagram 
in the [ijU-T jU plane, where a series of adjacent Mott- 
insulator lobes appears. Outside these lobes the cost of 
changing the total number of particles in the system by a 
finite amount vanishes in the thermodynamic limit. Con- 
versely, within the Mott lobes, adding or subtracting a 
single particle always costs a finite amount of energy, and 
the total boson population is commensurate to the size 
of the structure, N = uqM, where uq £ N. Hence, the 
left and right boundary of the filling-no Mott lobe are 
customarily evaluated as the ground-state energy gaps 
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where ii'Mott (j^o , t) , i?hoic("-o, t) and impart (»^o , t) are the 
ground state energies of H when the total population is 
npM noM — 1 and no A/ + 1, respectively. Following 
Ref. l34l in the sequel we term the latter two defect states 
and label "def" quantities relevant to either of them. 

In Ref. |3J the quantities in Eqs. Q are quite satis- 
factorily approximated by a third-order strong coupling 
perturbative expansion, the perturbative parameter be- 
ing r. Such results carry an explicit dependence on the 
hopping matrix t and, quite interestingly, on its maximal 
eigenvalue and eigenvector. This means that the influ- 
ence of topology, algebraically described by t, can be in 
principle singled out explicitly in the perturbative ex- 
pansion of ^L and Mr.- However, the perturbative results 
reported in Ref. JS^ apply only to homogeneous bipartite 
lattices, whose hopping matrix is such that X]/i=i ^jh = z 
and {t'^^^^)hh — 0. Actually, since such expansions are 
limited to the third order, the latter constraint can be 
weakened, and only regular lattices featuring triangular 
loops have to be excluded, {t^)hh = 0. In this respect we 
mention that third order results have been reported for 
the particular case of the regular triangular lattice in the 
presence of a magnetic field i^ 



III. PERTURBATIVE RESULTS 

In the following we provide the complete third-order 
strong coupling perturbative expansions for the quanti- 
ties in Eqs. (0), obtained including those terms arising 
from the presence of triangular loops and topological in- 
homogeneities in the discrete structure underlying the 
BH model. Such terms, pertaining to the third pertur- 
bative order, are denoted "<]" and "x", respectively. Quite 
interestingly, the latter depends not only on the maximal 
eigenpair of the hopping matrix, but on its entire spec- 
trum tik = Afeffc, Ai = maxfe(Afe). 

For commensurate fillings, N = n^M, the energy of 
the ground state is 
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where E^^^^^ 
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no(no - 1) - ^no(no + l)tr(i2) is 
the previously reported result9A, while f^ott = u^ ('^^ "'" 
l)(2no -I- l)tr(i'^) takes into account triangular loops, but 
it does not affect the boundary of the Mott lobe, since 
an identical term appears in the energy of defect states. 
Indeed, the latter has the form 

£;dcf = iJ^ott + ^def + r' (f Sou + Cf + ^dcf ) + o{r^) (4) 
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is the correct third-order result when triangular loops 
and hopping inhomogeneities are absent^ In detail, in- 
troducing the matrices CjTi = <5j7i S^^i i?,, (i'^')j7i 
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TABLE I: Values of the coefficients in Eqs. ®-@. Their 
correctness has been verified on small arrays (M < 10) by 
checking Eq. Q against exact numerical diagonalization. 
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where def denotes the defect state complementary to def. 
As it is clear from tabled the previously unreported^ 
third-order terms, 
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appearing whenever triangular loops are present, and 
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arising from topological inhomogeneity, can contribute 
significantly to the perturbative expansion of -Bdcf, and 
hence of Eqs. |(2Jl. Note that, similar to the previously 
known terms, ^"^ E^^f depends on the maximal eigenvector 
of t only. Conversely S^^f depends on all of the eigenvec- 
tors ffc and eigenvalues Afc of t. Of course £^^f vanishes on 



0. 
such 



structures lacking of triangular loops, where {t 
Likewise, S^^f vanishes on structures where C — zl, 
as regular lattices}^ since f /j • f ^ = 5hk- 

In the following we sketch the key steps involved in 
carrying out the results in Eqs. Q-®- As we mention 
above, the ground state energies appearing in Eqs. (0) 
are obtained restricting the BH Hamiltonian to the rele- 
vant number eigenspace and regarding the hopping term 
V = — X]j7i o.jtjhaf^ as a perturbation on the atomic un- 
perturbed Hamiltonian, Hq = 'tX^i^'jC^ ~ !)■ Ac- 
cording to standard perturbation theory, the ground 
state oi H = Hq + tV and the relevant energy are 
expanded in powers of the perturbative parameter r, 
*) = I]„'^"IV'™), E = Yjn'^'^^n- The perturbative cor- 
rections, obeying the hierarchy of equations 
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are then expanded on the eigenfunctions of iJo- Note 
that the ground state of the latter is non-degenerate 
for commensurate boson populations, N — uqAI, and 
Af-fold degenerate in the number eigenspaces relevant 
to defect states, N = uq M ± 1. Hence we conve- 
niently set IV'n) = Z^f^i^*" 1^) + Z^qT^"''!")' where 
the quantities denoted by Latin letters refer to the D- 
dimensional lowest eigenspace of Hq, while those denoted 
by Greek letters refer to the remaining eigenstates. That 
is Dmou = 1, -Ddof = M and Ho\x) = Ux\x), with 
Uk — u < Ua Va, k. Furthermore, for defect states, 

i(l)l-l/2 Y^M ^_ ^dcfl 
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we set \k) = [C)^J,Y 

the j-th component of f^, 5^°''' = 

l"o)) ~ ®j=i l"o)j is the non-degenerate ground state 
of Hq in the commensurate- filling case, TV = Muq. This 
makes the first perturbative corrections quite straight- 
forward. Indeed, one gets 1-00) = |l)i eo = "" and 
ef"^ — — C^JAi, where "1" henceforth stands for fc = 1. 
The subsequent corrections are obtained by projecting 
Eq. H1U|) on suitably chosen eigenvectors of Hq, with the 
assumption that || X;™^i T"|V'„)f = 1 + o{t"'+^). 

Here we focus on the third order correction for de- 
fect states, containing the previously unreported terms 
in Eqs. I© and I®. Projecting Eq. ifTUIl with n = 2 
onto l^o) = |1) one gets e' 



r.l"^-(a|y-eff|0i)/(u 



dcf 
3 



E„(l|y|«)rr, where 

2 — \ui\v — ci iv^i// v" ~ ""q) is in turn obtained set- 
ting n ~ I and projecting onto \a). Finally, projecting 
Eq. (|10|l with n — and n = 1 onto \a) and \k), respec- 
tively, gives the coefficients for the first order correction, 
F^"^ = {a\V\l)/{u-Ua) and 
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Putting these results together gives 
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whose first term provides the third order correction in 
Eq. ©, as well as those in Eq. Q. The second term 
in Eq. (|12|l is instead responsible for the "topological" 
correction in Eq. I®. Note that the latter is the lowest 
order term where the degeneracy of the ground state of 
Hq plays an explicit role. Indeed, the |fc)'s with fc > 1 do 
not contribute energy corrections at lower orders. In the 
non-degenerate case, N = Mno, the topological correc- 
tion is absent also at the third order, the counterpart of 
Eq. I|12|) featuring the first term only. 

In the case of a physically significant quantity such as 
the local density of states, the "topological terms" G\ ' 
in Eq. (|ll|l provide a perturbative correction at an order 
as low as the first. Indeed, a straightforward calculation 
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FIG. 1: (color online) A: Four fold star network; B: first Mott 
lobe for the four-fold star network (dark) and the ID lattice 
(lighter hue, for comparison); Main plot: Local density for 
the ground state of the four-fold star network as a function 
of the chemical distance from the central site, k. Dark and 
light hues refer to particle and hole defect states, respectively. 
Solid and dashed lines refer to the zeroth and first order per- 
turbative results (r/[/ = 0.1). 
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FIG. 2; (color online) Three different 2D structures and, in 
the same hue, the first Mott lobe thereof. For the inhomoge- 
neous Kagome lattice (C.) we choose the hopping matrix so 
that the average hopping amplitude is one. In detail tuh is 
either 0.8 (long links) or 1.2 (short links). 



IV. CONCLUSIONS 
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where (Tpart = — fhoie = 1- A comparison with ex- 
act numerical diagonalization on small arrays (Af < 10) 
shows that Eq. ((T!^ provides satisfactory results. A sim- 
ilar correction appears also in the two-point correlation, 
(^laj-a^lV'), which we do not discuss here. 

The first order correction to the local density aris- 
ing from topological inhomogeneity is clearly shown in 
Fig. n] where we considered the case of a four-fold star 

38.39 



networki&iS^ (see figure caption for details). 

In Figure 121 we compare the first Mott lobes for three 
different 2D structures, namely the Euclidean lattice 
(A.), the homogeneous (B.) and an inhomogcneous (C.) 
Kagome lattice^^. Note that the first two (homogeneous) 
structures have the same coordination number, z = 4, 
so that the only difference between their Mott lobes, 
Eqs. ©J comes from the presence of triangular loops in 
the latter. The inhomogcneous Kagome lattice has the 
same coordination as A. and B., but it features two kinds 
of hopping amplitudes and three kinds of sites (see figure 
caption for details). Hence its Mott lobes get contribu- 
tions from both Eq. © and Eq. (|5J|. 



In summary, we report the complete third order ex- 
pansion for the ground state of the BH model on a com- 
pletely arbitrary structure. We show that in the gen- 
eral case two kinds of previously unreported'^^ terms ap- 
pear. The first kind arises from the presence of triangular 
loops, and contributes to the ground state energy of both 
commensurate- filling and defect states. Due to exact can- 
celation only the latter influences the boundaries of the 
Mott lobes. The second kind of contribution occurs in the 
defect states of systems characterized by inhomogcneous 
hopping. As we mention above, these include superlat- 
tices with constant coordination or hopping amplitude, 
such as those considered in Fig. |2](C.) or in Ref. 18, re- 
spectively. Quite interestingly, hopping inhomogeneity 
produces perturbative terms that depend on the entire 
spectrum of the hopping matrix rather than on its max- 
imal eigenpair. Furthermore, other than a third order 
correction to the ground state energy, it gives rise to a 
first order correction in the local density of bosons. We 
remark that, following the procedure sketched in Ref. l35l 
our results can be adapted to field-tuned transitions in 
a superlattice of Josephson junctionsi^^ More in general 
they supply a basis for a comprehensive study aimed at a 
deeper understanding of the role of topology in quantum 
phase transitionsi^. 
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